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Utilizing the three-fireball picture within the quark combination model, we study system- 
atically the charged particle pseudorapidity distributions in both Au+Au and Cu+Cu 
collision systems as a function of collision centrality and energy, -^/sjvjv = 19.6, 62.4, 130 
and 200 GeV, in full pseudorapidity range. We find that: (i)the contribution from leading 
particles to dN ch /dn distributions increases with the decrease of the collision centrality 
and energy respectively; (ii)the number of the leading particles is almost independent of 
the collision energy, but it does depend on the nucleon participants N par t \ (iii)if Cu+Cu 
and Au+Au collisions at the same collision energy are selected to have the same N par t, 
the resulting of charged particle dN/drj distributions are nearly identical, both in the 
mid-rapidity particle density and the width of the distribution. This is true for both 
62.4 GeV and 200 GeV data, (iv)thc limiting fragmentation phenomenon is reproduced, 
(iiv) we predict the total multiplicity and pseudorapidity distribution for the charged 
particles in Pb+Pb collisions at y^sjvjv = 5.5 TeV. Finally, we give a qualitative analysis 
of the N c h/ < Npartfi > and dN^/dr// < N par t/2 > |^a;o as function of ^/sjvjv an d 
N par t from RHIC to LHC. 
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1. Introduction 

In relativistic heavy ion collisions at RHIC energies, the charged particles are pro- 
duced copiously in vacuum. The number of charged particles per unit pseudo- 
rapidity dN c h/drj 1 and in particular its dependence on some variables, such as ra- 
pidity, collision centrality and energy, are the important observables, from which 
a lot of information about the hot and dense matter created in collisions can be 
extractecfl^l From the pseudo-rapidity density and the transverse energy per par- 
ticle, one can determine via Bjorken method the initial energy density of the fireball 
which can provide one piece of evidence for the deconfmement phase transition. In 
the fragmentation region, the charged particle production, in general, is thought to 
be distinct from that at mid-rapidity, although there is no obvious evidence for two 
separate regions at any of the RHIC energies. The pseudo-rapidity density dN c h/dr) 
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in forward rapidity region carries some information of leading particles produced 
in collisions The experimental data about the charged-particle pseudo-rapidity 
density in both Au+Au and Cu+Cu collision systems have been presented by the 
PHOBOS collaboration QMS the PHENIX collaboration 14 , and the BRAHMS 
Collaboration HSllSI. The data for the scaled and shifted pseudo-rapidity distribu- 
tion dNch/dr)' /(Np ar t/2), exhibit the limiting fragmentation phen omeno n in both 
Au+Au and Cu+Cu collisions at different energies and centr alitie J 1 ^ | 1 8 | 

Recombination of partons HSHSI ; Partonic coalescence E3H2H an( j QCM ESHSI] 
have been made to described many observations. In our previous work using 
a Gaussian-like shape rapidity dist ributio n for constituent quarks as a result of 
the Landau hydrodynamic evolutiorPSEIS we have presented the pseudo-rapidity 
distributions of charged particles in Au+Au collisions as a function of collision 
centrality and energy. The calculation results are in good agrement with the data 
in central collisions. In peripheral collisions, our predictions are slightly lower than 
data in high rapidity range. The reason may be that we have not considered the 
contribution of leading particles. In present work, taking into account the leading 
particle influence, we apply a three-fireball picture B l 34 | 35 | ^ Q describe the evolution 
of the hot and dense quark matter produced in collisions, and obtain the rapidity 
distribution of the constituent quarks just before hadronization. Then let these 
constituent quarks combine into initial hadrons according to a quark combination 
rule, and allow the resonances in the initial hadrons to further decay into final 
hadrons with the help of the event generator PYTHIA 6.1 ESI. 



2. Three-fireball picture and quark combination model 

In this section, we introduce the three-fireball picture, which is used to describe the 
rapidity distribution of quark and antiquarks just before hadronization. In addition, 
we briefly introduce the QCM which describes the hadronization of these quarks 
and antiquarks produced in collisions. 



2.1. Three-fireball picture 

It is known that the nucleus-nucleus collisions at RHIC energies are neither fully 
stopped nor fully penetrated. As the incident nuclei penetrate through the target 
nuclei, the most of the collision energy is deposited in collision region to form a 
big central fireball, and the penetrating quark matter forms two small fireballs, 
i.e. target and projectile fireballs, in forward rapidity region. The charged hadron 
pseudorapidity distribution is the total contributions from the three fireballs. 

The big central fireball which contains the main part of collision energy, controls 
the rough shape of charged particle pseudorapidity distribution (width and height). 
Relativistic hydrodynamics has successfully described the evolution of system before 
hadronization. Here we use a Gaussian-type rapidity dis tributio n for constituent 
quarks as a result of the Landau hydrodynamic evolution HUHUS 
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/M-^-Pt-^). (i) 



where 



a 



2 2cj ,E^Js NN 



H^r^)- (2) 



1 - c* v 2 m 



Here, E is the effective energy offered by per participant pair, and it is used to 
produce the central fireball. m p is proton mass and c s is the sound velocity. e c is the 
energy in the volume of a free hadron at hadronization. All quarks and anti-quarks 
in the central fireball are within the rapidity range [—y m ax,ymax], 



c s , E^Js NN 
l + ci 2 m v e c 



The average constituent quark num ber in the big central fireball can be obtained 
from a simple quark production model E913QI 

(N q ) = 2[(a 2 + PE) 112 - a](AWt/2>, (4) 

where the parameter (3 ~ 3.6 GeV, and the parameter a — /3m — i, m is averaged 
quark mass and it is taken to be 0.36 GeV, they are the same with Ref. [29] 

The two penetrating fireballs mainly consist of leading light quarks. We also 
adopt a Gaussian type rapidity distribution for leading quarks 

r(v) i gW) «*p(-^) (5) 

J[y> (N q{T /P)) dy ^? ' {0> 

where N q (r/p) is the total quark number in the penetrating target and projectile 
fireballs, yo — ± yt """ 7 ' 2 f i ''"° x is the rapidity center for target and projectile fireballs 
respectively. Rapidity distribution range of quarks in the center of mass frame is 
y € [-Vbeam, -Vmax] for target fireball and y e [y m ax, Vbeam] for projectile fireball, 
respectively. In this work, the spectrum width of penetrating fireballs is taken to be 
a' = 0.18. 

The total energy of the three fireballs in nucleus-nucleus collisions is 

y/SNN (Np ar t/2), 

E ( t+p) = (Vijw- E)(N part /2), (6) 

where E( T+ p^ is total energy of the two penetrating fireballs. The average number 
of quarks in penetrating projectile and target fireballs (N q ^ T+P ^) is determined by 
E{T+py 

<W = 5§f. P) 
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(E q ) is the average energy of each quark in the penetrating projective/target fire- 
balls, and it can be written as: 

(E q )= J m T cosh(y)f'(y)dy, (8) 

Umax 

where rriT — \J m 2 + p\ is the transverse mass of leading quarks. The transverse 
momentum pt of leading quarks is approximately taken to be 0.25 GeV, one third 



of the value of net-proton at forward rapidity y r- ^ ' 



3. The quark combination model 

The QCM was first proposed for high energy e + e~ and pp collisions and recently it 
was extended to ultra-relativistic heavy ion collisions EH2l The model describes the 
production of initially produced ground state mesons (36— plets) and baryons (56 — 
plets). In principle the model can also be applied to the production of excited states 
™l These hadrons through combination of constituent quarks are then allowed to 
decay into the final state hadrons. We take into account the decay contributions of 
all resonances of 56 — plets baryons and 36 — plets mesons, and cover all available 
decay channels by using the decay program of PYTHIA 6.1 The main idea is 
to line up N q quarks and anti-quarks in a one-dimensional order in phase space, 
e.g. in rapidity, and let them combine into initial hadrons one by one following a 
combination rule (see section 2 of Ref. 32 for a short description of such a rule). We 
note that it is very straightforward to define the combination in one dimensional 
phase space, but it is highly complicated to do it in two or three dimensional phase 
space 1441. The flavor SU(3) symmetry with strang eness s uppression in the yields of 
initially produced hadrons is fulfilled in the model EDI43J 



4. Centrality, system size and energy dependence of 
charged-particle pseudo-rapidity distribution 

In this section, we will study the system size, energy and centrality dependence 
of pseudorapidity distribution of charged particles, and multiplicity distribution 
in mid and forward rapidity range respectively in relativistic heavy ion collisions. 
Moreover, if the particle production mechanisms in A+A collisions at RICH and 
LHC are the same, we predict the energy and centrality dependences of A 



<N part /2> 



nnH dN '= h I n 

anC1 drj<N part /2> \n~0- 

From Eq. (J6j) , we can see that there is only one free variable, i.e. E or E(t+p), 
which should be determined from the experimental data. In the present work, we 
determine the effective energy E for central fireball by fitting the pseudorapidity 
density d ^ h \ v ~q in Au+Au collisions at ^/s = 130 GeV. Then the E^ T+P ^ and 
leading quark number (N q rr+p)) can be naturally obtained. Using this method, we 
get the leading quark number in different collision centralities and parameterize it 



Centrality, system size and energy dependences of charged- particle pseudo-rapidity distribution 5 




Fig. 1. Charged hadron pseudo-rapidity distributions for different ccntralitics in Au+Au colli- 
sions at v / sjvjv = 200 GeV. The lines are our results and the points are data taken from PHOBOS. 

as the function of nucleon participants N part 

N q{T+ P } = -84.44 + 35.82 * N Q pa \ t . (9) 

At other energies, basing on the relation between the number of leading quarks 
and the centrality, we can get E, and the pseudorapidity distribution in different 
centralities within a quark combination model. 

Applying Eq. (J9j) to other RHIC energies, we firstly calculate the charged par- 
ticle pseudorapidity distributions for different central ities in Au+Au collisions at 
i/sWn = 200 GeV. The data are taken from PHOBOS EH2] The resu i ts are shown 
in Fig. [TJ The dashed lines are the contribution from the central fireballs. The dot- 
ted lines in the forward pseudorapidity range show the contributions of penetrating 
target and projectile fireballs respectively. The solid lines are the total contribution 
of the three fireballs. One can see that our results are in good agreement with the 
data. In addition, we find that the contribution from leading particles to dN^/di) 
distributions increases with the decrease of the collision centralities. We also give 
the results of dN c h/dr) distribution at v /saT/v=19.6, 62.4 GeV, and they are shown 
in Fig. [2] The data are taken from PHOBOS EEH The agreement of calculated 
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Fig. 2. Charged hadron pseudo-rapidity distributions for different centralities in Au+Au collisions 
at y/Sffpf =19. 6, 62 .4 GeV. The lines are our results and the points are data taken from PHOBOS. 



results with the data is also satisfactory. The Fig.[3]shows the charged-particle pseu- 
dorapidity distributions in most central Au+Au collisions at y / s]vjv=19.6, 62.4, 130 
and 200 GeV. The results indicate that the contribution from leading particles to 
dN c h/drj distributions increases with the decrease of collision energy. 

Recently, PHOBOS Collaboration have presented the data on charged-particle 

-jo 

pseudorapidity distributions in Cu+Cu collisions at ^/sjvat=62.4, 200 GeV 1 ^. The 
other goal of this paper is to investigate the systematic dependence of particle 
production in nuclear collision at RHIC energies, in terms of overall dN/dr/ dis- 
tributions. We apply Eq. ([9]) to Cu+Cu collisions, and give the charged-particle 
pseudorapidity distributions as a function of centrality at y / J/v r iv=62.4, 200 GeV. 
The results are shown in Fig. |4l and compared with the data. As we can see, our re- 
sults are roughly consistent with the experimental data, but are slightly lower than 
the data in the mid-rapidity range, especially for the central Cu+Cu collisions. The 
reason may be that we did not consider the difference between the collision ge- 
ometry in Cu+Cu and Au+Au, even at the same N part in Eq. ([5]), especially for 
the central Cu+Cu collisions. The two system have same shape, which indicate a 
similarity particle production mechanism between Au+Au and Cu+Cu, the results 
is the same as Ref. SSI Following that, we think the Au+Au and Pb+Pb collisions 
have the same particle production mechanism, so we predict the pseudorapidity dis- 
tribution and charged particle multiplicity in Pb+Pb collisions at y / J/v r /v=5.5 TeV 
in the followings. 

The Large Hadron Collider (LHC) is scheduled to begin operation in May 2008. 
The most pressing issue for the early days at the LHC is to establish the global fea- 
tures of heavy ion collisions. This involves the estimation of the inclusive charged- 
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particle yield and the charged pseudorapidity distribution and so on. In this work, 
extending Eq. ([9]) to LHC energy, we can predict the charged particle pseudora- 
pidity distribution as a function of centrality in Pb+Pb collisions. As an example, 
we calculate the most central Pb+Pb collisions at ^/snn — 5500 GeV in Fig. [5] 
The total charged-particle multiplicity is about 18170, and pseudorapidity density 
dN c h/dri |i^i<o.5 is about 1630. Note that the sound velocity is taken be 1/3. To 
separate the trivial kinematic broadening of the dN c h/drj distribution from the more 
interesting dynamics, we also study the scaled, shifted pseudorapidity distribution 
dN c h/drj' I (Np ar t/2} , where rf = rj — yteam , in most central Au+Au, Cu+Cu and 
Pb+Pb collisions at different energies. The calculation results are shown in Fig. [6] 
We observe that the data at various energies and systems fall on a common limiting 
curve. 

The centrality dependence of N c h/ < N part /2 > and dN c h/dr]/ < N part /2 > 
\ n ~o at y/sNN— 19.6, 62.4, 130, 200 GeV in Au+Au collisions (open symbol) com- 
pared with the data (solid symbol) taken from PHOBOsQ^ and the prediction at 
5500 GeV in Pb+Pb collisions (the solid line) are shown in Fig. [7l 

Fig. [8] shows the cm. energy dependence of N c h/ < N part /2 > and dN c h/drj/ < 
Npart/2 > \ v ;zo from RHIC to LHC predicted by the QCM (the lines). It shows that 
the N ch / < Np art /2 > and dN ch /dr\/ < N par t/2 > \„ao from RHIC to LHC can 
grow at logarithmically with y/sNN- 
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Fig. 4. Particles pseudo-rapidity distributions for different centralities in Cu+Cu collisions at 
v /sjvjv=62.4, 200 GeV. The lines are our results, and the data are taken from PHOBOS. 



5. summery 

Within a combination model, we study the charged particle pseudo-rapidity distri- 
butions in Au+Au and Cu+Cu collision systems as a function of collision centrality 
and energy (^/snn — 19.6,62.4,130 and 200 GeV), in full pseudo-rapidity range. 
We use a toy model, i.e. three fireballs, to describe the evolution of the hot and 
dense quark matter produced in collisions. The big central fireball which carries 
the main part of collision energy controls the rough shape of the charged particle 
pseudo-rapidity distribution. We apply the Landau relativistic hydrodynamic model 
to describe the the evolution of highly excited and possibly deconfined quark matter 
created in the big central fireball. As a result, we obtain a Gaussian- type rapidity 
spectra of constituent quarks before hadronization. The other two small fireballs 
in foreword rapidity carry the information of the leading particles. We also use a 
Gaussian-type rapidity spectra of constituent quarks before hadronization. Then 
we use our combination model to describe the hadronization of initially produced 
hadrons including resonances, whose decays are dealt with by the event generator 
PYTHIA 6.3^^1 Firstly, by studying the contribution of leading particles to charged- 
particle pseudo-rapidity distribution in Au + Au collisions for different centralities 
at 130 GeV, we extract the centrality dependence of the average number of leading 
quarks from the data. Then we extend it to other RHIC energies. We calculate the 
charged particle pseudo-rapidity distributions in both Au+Au and Cu+Cu collision 
systems as a function of collision centrality, at -^/s/vaF = 19.6, 62.4 and 200 GeV, 
in full pseudo-rapidity range. The calculation results are in good agreement with 
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Fig. 5. Charged particles pseudorapidity distributions calculated by QCM in most central Pb+Pb 
collisions at y/s = 5.5 TeV. 



data. To separate the trivial kinematic broadening of the distributions of the pseudo- 
rapidity density from more interesting dynamics, we compute the scaled and shifted 
pseudo-rapidity density distributions dN c h/drj' / (N part /2) with rf = 77 — ybeam at 
collision energies 19.6, 62.4, 130 and 200 GeV. The good agreement with data is 
found. Furthermore, we predict the total multiplicity and pseudo-rapidity distri- 
bution for the charged particles in most central Pb+Pb collisions at ^/snn = 5.5 
TeV. Through investigating detailed dN ch /drj distributions, we find that: (i)Thc 
contribution from leading particles to dN c h/drj distributions increases with the de- 
crease of the collision centrality and energy respectively; (ii)The number of leading 
particles is, independent of collision energy, only a function of nucleon participants 
Npartior the same system; (iii)If Cu+Cu and Au+Au collisions at the same colli- 
sion energy are selected to have the same N part , the resulting of charged particle 
dN/drj distributions are nearly identical, both in the mid-rapidity particle density 
and the width of the distribution. This is true for both 62.4 GeV and 200 GeV 
data. (iv)The limiting fragmentation phenomenon is reproduced. Furthermore, we 
predict the total multiplicity and pseudorapidity distribution for the charged parti- 
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Fig. 6. The scaled, shifted pseudo-rapidity density at v /s]vjv=19.6,62.4,130,200 GeV in most cen- 
tral Au+Au,Cu+Cu,and Pb+Pb collisions at v / 5=5.5TeV. The lines are our results, the symbols 
are data taken from PHOBOS. 

cles in Pb+Pb collisions at ^Js^n = 5.5 TeV, and find the N c h/ < N part /2 > and 
dN c h/dr)/ < Np ar t/2 > \ v ~ n from RHIC to LHC can grow at logarithmically with 

y/SNN- 
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